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We present the results of an experimental and theoretical study of the electronically excited 
(1) ^E^ state of *^Rb2 molecules. The vibrational energies are measured for deeply bound states 
from the bottom up to v' = 15 using laser spectroscopy of ultracold Rb2 Feshbach molecules. The 
spectrum of each vibrational state is dominated by a 47 GHz splitting into a Og and Ig component 
caused mainly by a strong second order spin-orbit interaction. Our spectroscopy fully resolves the 
rotational, hyperfine, and Zeeman structure of the spectrum. We are able to describe to first order 
this structure using a simplified effective Hamiltonian. 
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I. INTRODUCTION 

Progress in the field of ultracold atomic and molec- 
ular gases has always been strongly linked to develop- 
ments in molecular spectroscopy. Photoassociation spec- 
troscopy, for example, has been important for studies of 
ultracold atomic collisions and for production of ultra- 
cold molecules [H-Il]- In 2008, after carrying out spec- 
troscopic searches, several groups managed to produce 
cold and dense samples of deeply bound molecules in 
well-defined quantum states |5l4lO|. For this, a variety 
of clever optical transfer and filtering schemes were de- 
veloped which involved electronically excited molecular 
levels. These levels had to be properly chosen for high 
efficiency and selectivity of molecule production. This 
will also be the case for future experiments involving cold 
collisions [III , ultracold chemistry (l3l - [T5| , and test- 
ing fundamental laws via precision spectroscopy 
A detailed understanding of the excited molecular poten- 
tials is therefore necessary. Very recent work [l^ inves- 
tigated the spin-orbit-coupled and Ir'Hu states of 
Cs2. In other work a detailed analysis of weakly hound 
Rb2 levels of the excited Ig state close to the 5S'i/2+5Pi/2 
dissociation limit is currently under way [2^ (see also re- 
lated work in Fig. 13 of ij). 

Here we present measurements and analysis for deeply 
hound {v' = 0...15) levels of the (1)3^+ (5Pi/2 + 5S'i/2) 
state of ^^Rb2. This state is relevant for the production 
of deeply bound molecules in the a'^I]+ state via stimu- 
lated Raman adiabatic passage [1, Hlj. The a^S+ state 
is the energetically lowest triplet potential, giving rise to 
long lived molecules which are of interest for cold colli- 
sion experiments. The levels of the a'^E+ state have been 
mapped out and identified in detail in a recent publica- 
tion p^ . 

The (1) ■^1]+ state is not easily accessible in conven- 
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FIG. 1: Spectroscopy scheme. A tunable laser couples the 
Feshbach molecule level \i) to an excited level |«') with en- 
ergy E in the (1) "^E^ state. The excited state quickly de- 
cays spontaneously. Level positions in the (1) ^E^ potential 
are detected through resonantly enhanced loss of Feshbach 
molecules, h is Planck's constant and c is the speed of light. 



tional setups since the molecules in a Rb2 gas at ambient 
temperatures are found in the X^^'^ ground state. Elec- 
tric dipole transitions from this state to the (1) '^E+ state 
are forbidden by the g u symmetry selection rule. As 
a consequence it has only quite recently become possi- 
ble to explore in detail the Rb2 (1) ^S^J" state. Lozeille 
et al. [23] performed photoionization spectroscopy of ul- 
tracold Rb2 molecules produced by photoassociation in 
a magneto-optical trap to resolve the large 0" — Ig split- 
ting of the vibrational levels. Mudrich et al. |24i] used 
pump-probe photoionization spectroscopy of Rb2 formed 
on helium nanodroplets to measure the vibrational pro- 
gression of deeply bound levels. Our work goes beyond 
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FIG. 2: Energy diagram of the Zeeman structure for some 
weakly bound Rb2 molecular states below the \ fa ~ i,nif^ — 
1) + I/a = b,mf^ — 1) atomic asymptote. Dots mark the 
levels and magnetic fields used for the spectroscopy. For sim- 
plicity, throughout the paper we report excitation frequencies 
V with respect to the two-atom threshold at G. These are 
obtained by subtracting the known atomic Zeeman shift 
and the known molecular binding frequency Vm from the ex- 
citation frequency vl- The inset zooms into the region near 
the Feshbach resonance at 1007.4 G, marked with an arrow 
and "FR". 

this as we fully resolve the rotational, hyperfine, and Zee- 
man structure of the Rb2 levels in the (1) state. 

The starting point for our experiments is an ultracold 
ensemble of weakly bound Rb2 Feshbach molecules in a 
well defined quantum level, which has contributions from 
both the (1) '^E+ and X^T,'^ states. A scanning laser with 
sub- MHz short-term linewidth drives a one-photon tran- 
sition to levels in the (1) state (see Fig. 1). We 
obtain loss spectra for various magnetic fields from 
to 1000 G. The data are well described by an effective 
Hamiltonian which contains terms for molecular rotation 
as well as spin-spin, hyperfine, and Zeeman interactions. 

The article is organized as follows. Section [IT] of our 
article presents the experimental setup for the Rb2 spec- 
troscopy. In section lnll we show typical measured spectra 
and discuss their main features. Section lIVI presents the 
model that wc use to describe our data. Section |V] ex- 
plains the detailed level structure of our data based on 
the model. We conclude with a summary and a short 
outlook in section IVll 

II. EXPERIMENTAL SETUP 

Feshbach molecules in state \i) are irradiated by a pulse 
of light from a continuous- wave laser the wavelength of 
which is slowly scanned over the range 1000 to 1050 nm 



(Fig.[T]). The laser beam has an 1/e^ intensity waist ra- 
dius of 130 /im at the molecular sample. It is linearly 
polarized along the magnetic bias field B (parallel to 
gravity) and thus can only induce tt transitions. The 
light pulse typically lasts for 50 ms and has a rectangular 
temporal shape. When resonant with an \i) — \v') transi- 
tion, our laser induces losses in the Feshbach population 
due to excitation to \v') and subsequent fast (ns) decay 
to unobserved states. For each data point, a new ensem- 
ble of Feshbach molecules has to be prepared, a process 
which typically requires 28 s. 

The preparation of the Feshbach molecules is described 
in detail in [2^ [23| • In brief, following laser cooling and 
evaporative cooling, we trap an ultracold atomic cloud 
of 3x10^ *''Rb atoms close to quantum degeneracy in a 
3D optical lattice. A sizeable fraction (« 1/3) of the lat- 
tice sites are doubly occupied and the atoms are trapped 
in the lowest Bloch band. Afterwards, we ramp over a 
Feshbach resonance at a magnetic field of 1007.4 G (1 G 
= 10^4 T) (see Fig.H]) [Hj to obtain 3 x lO^* Feshbach 
molecules with no more than a single molecule per lat- 
tice site. The lattice depth for the Feshbach molecules 
is deep enough to prevent the molecules from colliding 
with each other, and we observe lifetimes of up to a few 
hundred ms. 

In order to investigate the Zeeman structure of the 
levels in the l'^5]+ state, we carry out the spectroscopy 
at various magnetic fields between and 1000 G. Thus, 
after we have produced the Feshbach molecules, the mag- 
netic field is ramped down from 1007.4 G to its cho- 
sen value. When doing so, one has to keep track of 
the quantum state \i) of the Feshbach molecules as well 
as their binding energy as they both change in gen- 
eral with magnetic field. Fig. [2] shows a number of rel- 
evant weakly bound molecular energy levels which ex- 
hibit many avoided crossings. When ramping over such 
a crossing, the molecules could potentially end up in two 
different quantum levels. This would be problematic for 
the spectroscopy and is avoided as described below. 

A large fraction of our measurements are carried out 
at 986.8 G with a binding energy of 22.7MHzx/i (see 
inset of Fig. [5]). The corresponding Feshbach state is 
best described in the atomic basis and correlates to 
\v = 40,(/, = IJt = 1)/ = 2,1 = 0,^^ = 2,Af = 2) 
at low magnetic fields. Here v is the vibrational quan- 
tum number of the a'^^t state, and fa,fb are the total 
angular momenta of atoms a and b. The sum of these an- 
gular momenta, f = fa + fb, and the rotational angular 
momentum of the atoms 1 couple to form the total angu- 
lar momentum F. At high magnetic field, F is no longer 
a good quantum number; only its projection M onto the 
z— axis remains good. As an estimate of this effect, we 
note, for example, that at 986.8 G the expectation values 
[2^ for fa and fb become about 1.5. 

For magnetic fields below 986.8 G we use a molecular 
level which correlates to \v — 36, {fa — 2, /(, = 2)/ = 
2, ; = 0, = 2, M = 2) at low magnetic fields. This is 
the diagonal line in Fig. 2 going from the point {B = 0, 
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E/h « —3.6 GHz) to the Feshbach resonance at thresh- 
old. It exhibits several avoided crossings with other 
molecular levels. When ramping down the magnetic field, 
these avoided crossings are crossed using a adiabatic ra- 
diofrequency transfer method ^27i] . In order to count the 
molecules remaining after the spectroscopy pulse, we re- 
trace our path back to the Feshbach resonance at 1007.4 
G where the molecules are dissociated via a reverse Fesh- 
bach magnetic field sweep and imaged as atoms using 
standard absorption imaging techniques [5^ . 

The spectroscopy laser is either a Ti: sapphire or a 
grating-stabilized diode laser, both of which can be 
locked to an optical cavity using the Pound-Drever-Hall 
scheme. The optical cavity is stabilized with respect to 
an atomic *^Rb line. While the unlocked lasers typically 
drift a few MHz within one experimental cycle, the cavity 
lock leads to a stability better than 1 MHz. Locking was 
only necessary for resolving a few weak lines. The laser 
frequency was determined using a commercial wavemeter 
(HighFinesse WS7) within seconds after the laser pulse. 
The wavemeter has an accuracy of 60 MHz after calibra- 
tion, which is done daily using an *^Rb line. Over the 
length of only a few experimental cycles (5 minutes) it 
typically drifts less than 10 MHz which sets a precision 
limit on relative line positions when working with an un- 
locked laser. 

We choose the |/a = I.^ti/^ = 1) + I/a = b,mf^ = 1) 
dissociation threshold at G as the energy reference {E = 
0). We thus substract from the measured laser excitation 
energy hv^, both the bound state energy of the Feshbach 
level hvrn as well as the Zeeman energy hva of the free 
atoms (Fig. [2]), which are both well known [22| . 



III. EXPERIMENTAL OBSERVATIONS 

In a first set of experiments, we have mapped out 
the vibrational ladder of the (1)'^S+ state from ti' = 
to v' = 15 at low resolution (Fig.|3^). We used the 
Ti:sapphire laser with an available power of a few hun- 
dred mW such that we only observed broad lines with 
a typical width of several GHz. The magnetic field 
was set to 986.8 G. The vibrational ground state of the 
(1) ■^1]+ state has an excitation frequency of 281.1 THz 
with respect to ji), corresponding to a laser wavelength 
of 1067nm. The vibrational splitting between the two 
lowest levels, v' — Q and v' = 1, is about 1.2 THz. The 
solid line in Fig. [3^ is a quadratic fit in v allowing for a 
slight anharmonicity in the vibrational ladder. The data 
lead to the same Morse potential as given by Mudrich 
et al. [23|. We find that each vibrational level is a dou- 
blet with a splitting of about 47 GHz [sl] (see Fig.[3]b). 
The large error bars of several GHz reflect the crudeness 
of this first measurement where we do not resolve the 
substructure of each of the doublet components. 
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FIG. 3: Vibrational ladder for the (1) '^S^ state, a) Excitation 
frequencies v of the vibrational states with quantum number 
v' . The solid line is a quadratic fit to the data, b) Zoom 
into the vibrational levels which are split into a Og and a Ig 
component. The spectra are measured at low resolution and 
high power. A frequency offset Vo(v') is subtracted from the 
excitation frequency u for each data point. The offset Voiy') 
is a quadratic fit to the Ig data. Note the different energy 
scales of a) and b) . 



A. Splitting of the vibrational levels into Og and Ig 
components 

The 47 GHz splitting of the vibrational levels clearly 
cannot be explained by the rotational, hyperfine, or Zee- 
man interactions because they would be too small. Es- 
timating the molecular hyperfine and Zeeman energies 
from those for ^''Rb atoms, we expect such contributions 
to be at most 14 GHz. It turns out that the large split- 
ting comes from a strong effective spin-spin coupling of 
the electrons. Although there is direct spin-spin coupling, 
the main contribution is second order spin-orbit coupling, 
which is resonantly enhanced by the nearby (l)^Hg state 
[2^ . Experimentally, these two contributions cannot be 
separated. In its microscopic form, the effective spin-spin 
interaction Hamiltonian reads 



-C 



3(si ■ ri2)(s2 • ri2) - (si • S2)rf 



(1) 



where C is a constant, Si^2 are the electron spin oper- 
ators, and ri2 is the relati ve p osition vector of the two 
electrons. It can be shown [Ulll] that for E states, Hss 
can be simplified to 



i/.. -2A[(n -8)2-8 V3]. 



(2) 



Here, A is effective molecular parameter for the spin-spin 
interaction which has to be determined for the studied 
vibrational level, n is the unit vector along the inter- 
nuclear axis, and 8 is the total electronic spin operator 
with 8 = Si -|- S2. Since the term 8^/3 merely results in 
an overall offset, it will be ignored. Thus the spin-spin 
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interaction takes the form 



H,, = 2A(n • S) 



(3) 



A strong Hss couples the electronic spin to the internu- 
clear axis, making its projection S = 0, 1 a good quantum 
number. Thus, for our (1) state fi^A + E^^Eis 
also a good quantum number. Here A and D, are the pro- 
jections of the total electronic orbital angular momentum 
L and the total angular momentum J = L + R+S on the 
internuclear axis, where R is the rotational angular mo- 
mentum of the nuclei. The eigenvalues of Hss are 2AE^. 
This means that the splitting between the = and 
n = 1 states is 2A. As will become clear later, the more 
deeply bound component of the observed doublet struc- 
ture has Ig character and the other one has 0~ character 
(Fig.|3]b), where we use Hund's case (c) notation 



B. Spectra of the Og 

The Og and Ig states have a rich substructure which 
we are able to resolve by lowering the power of the laser 
to about 0.1 mW. Figures 2] and [5] show loss spectra for 
Ig {v' = 0, 13) and Og (v' = 13), respectively. While 
the Ig manifold is spread out over 12 GHz, the 0" man- 
ifold is narrower (3 GHz) and has fewer lines. The 12 
observed lines of the Ig manifold in Fig. 2] show no ob- 
vious pattern. The structure is the result of an interplay 
of rotational, hyperfine, and Zeeman interactions. It is 
one of the main goals of this work to understand this 
structure and to identify the individual lines. 

A Ig spectrum typically consists of roughly 300 points 
where each point corresponds to one production and mea- 
surement cycle. The lines that we observe vary markedly 
in width. This indicates a strong variation of the laser- 
induced coupling between \v') and Spectra for differ- 
ent vibrational levels are similar (Fig.|4]). 

Compared to the typical step size of 40 MHz in Fig. 21 
the 12 MHz natural linewidth of the molecular levels 
(given by twice the atomic linewidth) is relatively small. 
Thus, it is possible that some weak lines are not always 
detected, especially when they are located on the shoul- 
der of a power-broadened line. We thus carried out a 
number of scans with various step sizes, laser powers, and 
pulse times, testing for consistency and checking theoret- 
ical predictions. For the Ig spectra we found 18 lines in 
total (Section El. 

Our Og" spectrum is considerably simpler than the Ig 
spectrum. There are 5 lines arranged in a doublet-like 
structure with a splitting of w 2.5 GHz (Fig.[51). This 
splitting is due to rotation of the molecule [40l|, 



(4) 



where Byi is the rotational constant and h is the re- 
duced Planck constant. The hyperfine and Zeeman con- 
tributions are much weaker here. As we will show in 
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FIG. 4: The Ig spectra of two different vibrational levels 
at 986.8 G. For each spectrum, we subtract a frequency off- 
set Vo from the excitation frequency v such that the most 
deeply bound line is situated at GHz. a) v' — 13, i^o = 
294626.4 GHz. b) v' = 0, Uo ^ 281068.2 GHz. In these two 
scans, not all lines that we eventually observed are resolved. 



section IIVI this can be explained due to the vanishing of 
the projections of total spin S and total orbital angular 
momentum L of the electrons on the internuclear axis. 
The rotational constant By' 



IS 



B„, 



2m\ 



(5) 



where |f') is the ket for the vibrational wave func- 
tion, r is the nuclear separation, and u is the reduced 
mass. From previous investigations, [23l . |30| we expect 
Byi^is/h = 412 MHz. Due to the weakness of the hyper- 
fine and Zeeman interactions, the angular momentum J 
is conserved. Apart from an offset the rotational energy 
is then determined from Eq. (|4]) to be 37 1 



Er, 



By, J(J+ 1) 



(6) 



We observe the lines J = and J = 2 separated by 
6By,/h w 2.5 GHz (see Fig. 5). The rotational level J = 
1 is not accessible because total parity (under inversion 
of all electron and nuclear coordinates) has to change 
in the optical transition. The parity for the Feshbach 
molecule is (—1)' +1 while the parity of the states in 
the 0" manifold of the (1) state is (-1)"^+^ The 
J — 2 feature has a substructure of four lines which is 
due to residual hyperfine and Zeeman interaction. This 
substructure will be discussed in detail in section fVl hand 
in hand with the analysis of our theoretical model. 



IV. EFFECTIVE HAMILTONIAN AND 
EVALUATION OF ITS PARAMETERS 

In the following we present the diatomic model Hamil- 
tonian which we use to describe the observed energy lev- 
els within the v' = 13 vibrational manifold. Spectra in 
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FIG. 5: Scans of 0" (u' = 13) lines in the (1) 
a magnetic field of 986.8 G. The dotted line corresponds to a 
laser power of 0. 1 mW and a pulse duration of 200 ms while the 
solid line corresponds to a laser power of 0.1 mW and a 1ms 
pulse duration. The frequency offset is i/o ~ 294671.0 GHz. 
The line at u — Uo — GHz corresponds to J = 0, while the 4 
lines located around 2.5 GHz have J — 2. Due to saturation 
these 4 lines are not resolved for the 200 ms pulse. 



other low lying vibrational manifolds are similar (Fig. 
4), and can be described by the same Hamiltonian with 
slightly adjusted parameters. We can thus ignore elec- 
tronic terms and vibrational motion. The Hamiltonian 
can then be written, 



H = H., 



hf 



(7) 



Hss = 2A(n • S)^ is the effective spin-spin operator 
given in Eq. (3) which, as previously discussed, leads 
to the large splitting into the 0" and Ig components. 
H^ot = ByiH^/h^ is the Hamiltonian for molecular rota- 
tion from Eq. (4). The terms for hypcrfinc interaction 
ffhf, Zeeman interaction Hz, and finally spin- rotation in- 
teraction Hsr are described in the following. Since the 
goal in this article is to get a first understanding of the 
experimentally observed spectra, we will in general sim- 
plify the interaction and only take into account terms of 
leading order. 



A. The hyperfine interaction 

A general ansatz for the hyperfine Hamiltonian in 
Hund's case a) and b) is of the form [s^l 



fl"hf = a A I • n -I- - -c) I • S -I- c (I • n) (S • n). (8) 

3 

Here, I is the operator for the total nuclear spin. The 
first term describes the interaction of the electronic or- 
bital angular momentum with the nuclear spin. However, 
since A = 0, this term will not contribute. The Fermi 
contact parameter is 6^^, while c is called the anisotropic 
hyperfine parameter. For S-states, we expect c ^ bp 
f|32|. page 196). These two parameters could in princi- 
ple be calculated ab-initio, but we have used them as free 



fit parameters. We note that the total nuclear spin I is a 
good quantum number for the given Hamiltonian. 



B. The Zeeman interaction 

Because we carry out measurements at high magnetic 
fields, the Zeeman interaction plays an important role. 
The main contribution to the Zeeman interaction comes 
from the electrons while contributions from the nuclear 
spins and molecular rotation (as well as second order ef- 
fects treated in f36| ) are much smaller and are neglected 
here. Furthermore, since A vanishes in the (1) state, 
the Zeeman interaction due to the total orbital angular 
momentum of the electrons, /isL • B/?i is also negligible 
to first order. (Here, /is is the Bohr magneton). The 
only remaining term is 



Hz = M-BffsS • B/h, 



(9) 



where gs is the electron g-factor. We note that there 

is no free fitting parameter in the Zeeman Hamiltonian 

for adjusting the model to the measured data. While the 

contribution from the Zeeman interaction will in general 

be large for the 1^ fines, it will be small for the Or lines 

^ . . . . ^ 

because here the spin projection onto the internuclear 

axis E = 0. 



C. The spin- rotation interaction 

In principle, we could also include a spin-rotation in- 
teraction through the effective Hamiltonian 



(10) 



where N = L + R and is the spin-rotation coupling 
constant. However, since "fyi is typically a small frac- 
tion of the rotational constant By' |37| , the spin- rotation 
interaction only represents an insignificant correction to 
the energy levels in the present system and is thus ne- 
glected. 



D. Fit procedure and evaluation of molecular 
parameters 

According to our model Hamiltonian in Eq. ([7]) there 
are four adjustable parameters: the rotational constant 
Byi, the spin-spin splitting parameter A, the Fermi- 
contact parameter bp, and the anisotropic hyperfine pa- 
rameter c. As discussed before in section IIIIBl the ro- 
tational constant By/ should be close to 412 MHz based 
on previous work, and in agreement with our analysis of 
the Og spectrum in Fig. 5. The spin-spin parameter A is 
determined by the splitting of the Ig and 0" manifolds 
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to be 2A = 47 GHz. This leaves only bp and c as com- 
pletely free parameters. We determined all parameters 
from fits of the model to the experimental data using a 
nonlinear Levenberg-Marquardt method. For the calcu- 
lations, the Hamiltonian in Eq. ([7]) is expressed in terms 
of matrix elements in a Hund's case (aa) basis, where the 
basis states are of the form 

\A,s,j:,h,i2,nj,,nj,,F,nF,M), (ii) 

where ^If = A -I- S -I- il/^ -I- flj^ is the projection of 
the total angular momentum F on the internuclear axis. 
Ii = I2 ~ i/l are the nuclear spins of the two nuclei 
and ri/j^ , r2/2 are their projections onto the internuclear 
axis. We note that the basis set PT|) is chosen larger than 
necessary for our purpose. It can be conveniently used 
to investigate also Hamiltonians where the total nuclear 
spin I (I = Ii -I- 12) might not be a good quantum num- 
ber. For the analytical expressions of the matrix elements 
as a function of the quantum numbers of Eq. we 
refer the reader to [39|. The Hamiltonian matrix is then 
numerically diagonalized to obtain the eigenvalues and 
eigenstates. Included in this calculation are all hyperfine 
states in the = 13) electronic state with total an- 

gular momentum of up to -F = 10. Experimentally, we 
have only observed states with F < 7. The parameters 
bp and c are determined in terms of combination bp + ^c 
and bp — ic, which correspond to contributions diagonal 
and off-diagonal in S and ilj = O/^ -I- flj^ , respectively. 
The diagonal term can be directly read off from Eq. ([SJ . 



= {bp - -c) 1]/E + cf7/E 

= {bp + -c) rtij:. 

Our final analysis gives By' = 412 MHz, 2A = 47 GHz 
and &F + |c = 832 MHz. We find, however, that bp - 
is not precisely determined in our analysis. Experimen- 
tally, the parameter bp — could, in principle, be best 
determined from the 0" spectrum because the contribu- 
tion to the hyperfine interaction is purely off-diagonal 
with respect to S. However, the energy levels depend 
only weakly on this parameter, and our measurements 
indicate a value between (200 — 1000) MHz. Higher pre- 
cision than that reached in our measurements is required 
in order to better estimate {bp — ^c). 

V. ASSIGNMENT OF LINES OF THE 0" AND 
Ig SPECTRA 

With the help of our theoretical model, we can now 
identify the individual lines of the 0" and Ig spectra 
and understand the physical origin of the substructure 
of these spectra. Since the Feshbach molecules have to- 
tal magnetic quantum number Af = 2, and since we use 
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FIG. 6: Zeeman diagram of the 0^ {v' = 13) state. An offset 
Vo = 294600 GHz is subtracted from the excitation frequency 
V. The lines are calculations. Dashed lines indicate levels 
with _F > 3 that cannot be excited in our experiment at B = 
due to the selection rule on AF — 0, ±1. Table |I] gives the 
quantum numbers at G and 1005.8 G. The inset is a zoom 
into the closely spaced lines around 1000 G. 



an optical tt transition, we only observe excited levels 
with M = 2 and thus F > 2. From our discussion so far, 
we have the following good quantum numbers for our 
levels in the (1) ^Ej state: M = 2, A = 0; E = 0, 1; O = 
0, 1; 5 = 1. Furthermore the total nuclear spin / is a good 
quantum number in our Hamiltonian and as we will see 
in the next section can take the values 7 = 1,3. At low 
magnetic field, F is also a good quantum number. With 
these and the additional quantum number J, which (as 
we will show) is good for the 0" states and becomes good 
for high-rotational Ig states, we will be able to identify 
all lines. 



A. Og Zeeman spectrum 

Figure [6] shows experimentally observed energy levels 
for various magnetic field strengths along with the cal- 
culations (dashed and solid lines). The Zeeman effect is 
quite weak as levels shift no more than a couple hundred 
MHz over a range of 1000 G. This is consistent with our 
previous assertion in section Illb that the Zeeman inter- 
action (cx S • B) is small because the expectation value 
of the electron spin vanishes in any projection direction 
for a state with E = 0. An identical argument holds 
for the hyperfine interaction in Eq. (8). The observed 
residual splitting of the lines in Fig. 6 is mainly due to 
second order contributions of the Zeeman and hyperfine 
interactions. The 0" spectrum is dominated by the ro- 
tational splitting ByiJ{J + 1). This overall trend is well 
described by the theory (Fig. [5]). A slight discrepancy 
can be observed in the rotational splitting which seems to 
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lie somewhat outside the experimental error bars. Note 
that at G, experimental data are less accurate than 
for the other magnetic fields because of technical stabil- 
ity issues when we transfer the Feshbach molecules from 
1007 G to low magnetic fields. Uncontrolled magnetic 
field drifts lead to fluctuations in the transfer efficiency 
across the avoided crossings from shot to shot, which 
made the measurements noisy. 

We now study in detail the quantum numbers of the 
observed levels. Based on the exchange symmetry of the 
nuclei and other fundamental symmetries of the (1) 
state, one can show that the J = and J = 2 levels 
have either total nuclear spin / = 1 or 3. J and I couple 
to form the total angular momentum F. For J = we 
obtain F — 1,3, of which only = / = 3 is observed (the 
single line at « 294671 GHz in Fig. 6). This is because 
we only detect levels that have a total magnetic quantum 
number M = 2. Coupling J = 2 with / = 1 (J = 3) 
results in states with F = 2, 3 (i^ = 2, 3, 4, 5), where we 
have omitted F = 1 states, since again they cannot be 
observed. These six levels are plotted as lines in Fig. 6 
and also listed in tableUtogether with measured data. As 
can be read off from the expectation values in the table, 
for low magnetic fields the levels in the Og component are 
indeed well described by the quantum numbers J, /, F. At 
1000 G, F is not good anymore. 

In Fig. [HI two of the calculated lines are dashed. They 
correspond to F = 4, 5 which are not observable at OG 
because of the selection rule AF = 0, ±1. For larger mag- 
netic fields, the Zeeman interaction mixes states with dif- 
ferent F and the two levels that correlate with F = 4, 5 
at G become observable. As the magnetic field is in- 
creased, some lines cross. Because / is exactly a good 
quantum number for our specific Hamiltonian H, cross- 
ings between levels with different / quantum numbers are 
not avoided. 

From Fig. 6 as well as from Table I, we see that ap- 
parently the two 1=1 lines are not observed. Closed 
coupled-channel calculations 1:22, 38] show that our Fes- 
hbach molecules only have a few percent 1=1 char- 
acter. The selection rule AJ = thus makes it clear 
that 1=1 lines will be suppressed in the spectrum. 
We also know that the Feshbach state is a superposi- 
tion of electronic singlet and triplet states. For I = 
one finds from symmetry arguments (inversion symme- 
try, rotational symmetry, exchange symmetry, refiection 
symmetry) that the triplet component has nuclear spin 
/ = 1,3 while the singlet component has / = 0, 2. Fur- 
ther, for I = 0, the triplet component of the Feshbach 
state in the Hund's case (c) basis has J = 1, fi = while 
the singlet component has J = 0. Only the electronic 
triplet part of the Feshbach state will contribute to the 
optical transition to the (1)'^E+ 0" state, fulfilling the 
standard Hund's case (c) = $1 = selection rules 
A J = ±1, AS* = 0, and AI = 0. 





B = OG 




B 


= 1005.8 G 




{1} 


(F) 


— Vo 

(GHz) 


VE — Vo 

(GHz) 


(/) 


(F) 


Vt - Vo 

(GHz) 


Ve — Vo 

(GHz) 


0.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 


3 
3 
1 
1 
3 
3 
3 


3.0 
5.0 
3.0 
2.0 
4.0 
3.0 
2.0 


71.44 
73.84 
73.84 
73.92 
73.94 
74.06 
74.17 


71.14 

n.o. 

n.o. 

n.o. 

n.o. 
73.64 
73.95 


3.0 
3.0 
3.0 
1.0 
3.0 
1.0 
3.0 


3.0 
4.2 
3.7 
2.9 
3.3 
2.1 
2.8 


71.42 
73.83 
73.95 
74.00 
74.13 
74.27 
74.42 


71.63 
73.83 
73.92 

n.o. 
74.13 

n.o. 
74.40 



TABLE I: Calculated and measured OJ (v' = 13) excita- 
tion frequencies of the levels shown in Fig. [6] at G and 
1005.8 G. The offset from the excitation frequency v is Vo = 
294600 GHz. The subscripts t and e denote theoretical and 
experimental values. The expectation values for quantum 
numbers J, I and F are also given. All levels have S = 1, 
S = 0, and M = 2. Levels with 1 = 1 are not observed 
(n.o.). The absolute accuracy for the measured line positions 
at B = and 1005.8 G is about 200 MHz and 70 MHz, respec- 
tively. The accuracy for the relative positions of the lines, 
however, is in general higher, e.g. ~20MHz for the 3 = 2 
lines at B = 1005.8 G. 



B. Ig Zeeman spectrum 

The experimental Ig data at various magnetic fields 
together with the calculated energies are shown in Fig.[71 
There is good agreement of the overall structure of the 
calculated levels and the observed data. Our calcula- 
tions (Table [lT| indicate that the angular momentum J 
becomes a good quantum number for J > 3. Indeed, for 
J > 3, rotational lines belonging to different J are ener- 
getically well separated and the rotational splitting starts 
to dominate the structure of the observed spectrum. 

This can be understood as follows: A fast rotation 
of the molecular axis averages out the direction of the 
electron spin in the lab frame, because the electron spin 
is tightly coupled to the molecular axis (S = 1). This 
leads to an effective decrease of the Zeeman interaction. 
Similarly, a fast rotation will prevent the nuclear spin 
/ from coupling to the molecular axis (and therefore to 
the electron spin), as it cannot follow fast enough in an 
adiabatic way. In fact, calculated expectation values for 
ri/ (see Table ini are close to (typically |(0/)| < 0.5 
for J > 3) indicating a strong averaging out. For slow 
rotations (J < 2), in contrast, the hyperfine and Zeeman 
interactions are of the same order as the rotational split- 
ting. This leads to a relatively strong mixing between 
levels with J = 1 and J = 2. Note that for the Ig man- 
ifold only levels with J > 1 exist since the projection 
n = l. 

The energetically lowest state has an expectation value 
(|S1/|) = 2.6, which is close to its "ideal" value of 3 if the 
hyperfine interaction (cx Jl/S) was dominant over the 
rotational and Zeeman interactions. We can now ask how 
large this hyperfine interaction is in the limit of smallest 
rotation, e.g. J = 1. We can roughly estimate this by 
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FIG. 7: Zeeman diagram of the Ig {v' = 13) state. An offset Uo — 294600 GHz is subtracted from the excitation frequency v. 
The hues are calculations. Dashed lines indicate levels with _F > 3 that can not be observed in our experiment at S = due to 
the selection rule F = 0, ±1. The arrow points to an experimentally observed line that is not well described by the calculations. 
Table [III gives the corresponding quantum numbers at OG. 



looking at Fig. 7. At S = the rotational splitting 
of the J = 1 and J = 4 levels has to be 18 B^r — 7.2 
GHz. The observed splitting between the bary center of 
the J = 4 levels and the lowest line at v — ~ 26.34 
GHz is larger by about 2 GHz. This indicates that the 
hyperfine structure for J = 1 must be spread out over a 
range of roughly 4 GHz. 

We now investigate the number of lines that can in 
principle be observed. For J = 4 and / = 1,3 we can 
form 9 states with F > 2. For J = 3 we expect 8 states 
and for J = 1, 2 together we expect 10 states. These 
are the states which are listed in Table II. In contrast 
to the spectrum of 0~, we can observe levels with even 
as well as odd J, due to a two-fold degeneracy of each 
rotational level. One of the two levels has positive parity 
while the other one (that is observed in the experiment) 
has negative parity. 

The fact that we observe levels with a rotation up to 
J = 4 might at first be surprising, as we start from J = 1 
— in the Feshbach state and the selection rule is 
AJ = 0, ±1. However, this simply shows that J in the 
excited state is not a very good quantum number yet and 
the J w 4 levels have a J = 2 contribution. 

At zero magnetic field, the quantum number F for the 
total angular momentum is good. At this field we can- 
not observe levels with F > 3, due to the selection rule 
AF ~ 0, ±1. As in Fig. 6, these levels are drawn with 
dashed lines. Again, at larger magnetic fields levels with 
different F mix, and as a consequence more levels can be 
reached. The curves in Fig. [7] also display level crossings. 
Crossings are in general avoided as long as the levels have 



the same / quantum number. For our specific Hamilto- 
nian H, where / is exactly good quantum number, levels 
with different / do not mix with each other. 

The overall structure of the observed spectroscopic 
lines is well described with our model, which essentially 
only has a single free parameter, bp + |c. However, we 
find significant deviations of up to a few 100 MHz (e.g. 
note the data point next to the small horizontal arrow 
in Fig. 7). Such deviations clearly lie outside our exper- 
imental uncertainty, especially since for these cases rel- 
ative positions of neighboring lines are determined with 
a precision better than typically 30 MHz. On this level 
of accuracy, in general, we do not have good agreement 
with the theory. 

In order to achieve better agreement, one must include 
terms in the theory that we have so far neglected. For 
example, the Zeeman term oc L • B can contribute in sec- 
ond order. It could also contribute to first order if the 
projection A is not completely vanishing due to mixing- 
in of a A = 1 component from nearby states, e.g. via 
second order spin-orbit interaction. However, such inves- 
tigations would also require more detailed measurements 
and must be left for future work. 



VI. SUMMARY AND OUTLOOK 

In this article we carry out high resolution molecu- 
lar spectroscopy starting with an ultracold ensemble of 
^^Rb2 molecules. We are able to resolve the vibrational, 
rotational, hyperfine and Zeeman structure of deeply 
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= OG 


{J) 




1 /O \ 1 

\{ih)\ 


(F) 


f -i^o (GHz) 


1.2 


3.0 


2.6 


2.0 


26.34 


1.4 


3.0 


2.0 


3.0 


27.06 


2.7 


3.0 


1.2 


4.0 


28.36 


1.2 


1.0 


0.16 


2.0 


28.53 


2.0 


3.0 


1.0 


2.0 


29.00 


2.0 


3.0 


0.41 


3.0 


29.68 


1.9 


1.0 


0.27 


2.0 


29.98 


2.1 


3.0 


0.75 


5.0 


30.12 


2.3 


1.0 


0.031 


3.0 


30.13 


2.0 


3.0 


0.57 


4.0 


30.51 


3.0 


3.0 


0.75 


2.0 


31.77 


3.0 


3.0 


0.21 


3.0 


32.02 


3.0 


1.0 


0.11 


2.0 


32.31 


2.9 


3.0 


0.60 


4.0 


32.44 


3.0 


1.0 


0.025 


3.0 


32.45 


3.1 


3.0 


0.44 


6.0 


32.56 


3.2 


1.0 


0.069 


4.0 


32.57 


3.0 


3.0 


0.37 


5.0 


32.64 


3.9 


3.0 


0.13 


2.0 


35.04 


3.9 


3.0 


0.21 


3.0 


35.41 


4.0 


1.0 


0.014 


3.0 


35.54 


3.9 


3.0 


0.48 


4.0 


35.78 


4.0 


1.0 


0.13 


4.0 


35.81 


4.0 


1.0 


0.011 


5.0 


35.99 


3.9 


3.0 


0.54 


5.0 


36.07 


4.1 


3.0 


0.25 


7.0 


36.19 


4.0 


3.0 


0.30 


6.0 


36.22 



TABLE II: Calculated Ig («' = 13) excitation frequencies 
of the levels shown in Fig. [7] at OG. The offset from the 
excitation frequency v is Vo = 294600 GHz. The expectation 
values for quantum numbers J, I and F are also given. All 
levels have S = 1, E = 1, and AI — 2. Levels with I — 1 are 
not observed. 



Hamiltonian. In brief, the level structure for the 0" line 
is mainly determined by the vanishing spin component 
E = 0, which leads to a very small hyperfine and Zeeman 
structure and a good quantum number J. In contrast, 
for the Ig line, the hyperfine and Zeeman interactions 
are large for small rotations, but then are averaged out 
at larger rotations such that the rotational splitting ac- 
cording to J( J + 1) again determines the spectrum. De- 
spite the overall understanding of the level structure, we 
still observe systematic deviations between experiment 
and theory on the order of a few 100 MHz, which should 
disappear in a more refined model. We find from our 
data that the anisotropic part of the hyperfine interac- 
tion could essentially not be determined. In order to do 
this, the experimental data of the 0" line will have to be 
measured with higher precision in the future. 

It would be interesting to see whether the determined 
fit parameters for A (effective spin-spin interaction) and 
for the hyperfine contact interaction can be deduced from 
ab-initio calculations and the known atomic properties. 
Besides gaining a better insight into the level structure of 
the so-far relatively unexplored (1)'^S+ state, this work 
will be helpful for future cold molecule experiments where 
the deeply bound molecules need to be prepared with 
high efficiency in various well defined quantum states of 
the a'^S+ state. Levels in the (1) '^E^ state may then 
serve as an intermediate state-selective step in a two- 
photon optical transfer scheme from Feshbach molecules 
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bound states {v' < 15) of the excited (1) state. The 
accuracy of the measured lines is about 60 MHz, how- 
ever, the relative position of the lines with respect to 
each other is often determined much better and its ac- 
curacy reaches 20 MHz. The dominating feature of the 
observed spectra is the splitting of the vibrational levels 
into a Ig and a 0" component which can be understood 
as a strong effective spin-spin coupling of the electrons. 
We obtain a good understanding of the level structure to 
first order with a relatively simple effective model that 
only takes into account the most important terms of the 
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